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Kazal-type serine proteinase inhibitors 
in the midgut of Phlebotomus papatasi 

Leah Theresa Sigle, Marcelo Ramalho-OrtigaoA 

Department of Entomology, Kansas State University, Manhattan, KS, USA 

Sandflies (Diptera: Psychodidae) are important disease vectors of parasites of the genus Leishmania, as well as 
bacteria and viruses. Following studies of the midgut transcriptome o/Phlebotomus papatasi, the principal vector of 
Leishmania major, two non-classical Kazal-type serine proteinase inhibitors were identified (PpKzll and PpKzl2,). 
Analyses of expression profiles indicated that PpKzll and PpKzl2 transcripts are both regulated by blood-feeding in 
the midgut ofV. papatasi and are also expressed in males, larva and pupa. We expressed a recombinant PpKzl2 in a 
mammalian expression system (CHO-S free style cells) that was applied to in vitro studies to assess serine proteinase 
inhibition. Recombinant PpKzl2 inhibited a-chymotrypsin to 9.4% residual activity and also inhibited a-thrombin 
and trypsin to 33.5% and 63.9% residual activity, suggesting that native PpKzl2 is an active serine proteinase inhibi- 
tor and likely involved in regulating digestive enzymes in the midgut. Early stages q/Xeishmania are susceptible to 
killing by digestive proteinases in the sandfly midgut. Thus, characterising serine proteinase inhibitors may provide 
new targets and strategies to prevent transmission q/Xeishmania. 
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In arthropods, serine proteinases are involved in 
digestion, coagulation, phenoloxidase activation and 
other immune responses. Regulation of these enzymes 
by serine proteinase inhibitors is critical for maintain- 
ing homeostasis (Kanost 1999, Jiang & Kanost 2000, 
di Cera 2009). Several serine proteinase inhibitors have 
been identified in blood-feeding arthropods and linked 
to inhibition of thrombin and other components of the 
coagulation cascade to facilitate fluidity in the mouth 
parts and midgut following blood-feeding on a host [re- 
viewed by Tanaka-Azevedo et al. (2010)]. Many of these 
thrombin inhibitors belong to the family of Kazal-type 
serine proteinase inhibitors. 

The first Kazal-type thrombin inhibitor identified in 
a haematophagous insect was from Rhodnius prolixus 
(Friedrich et al. 1993). Since then, proteins containing 
Kazal-type domains have been identified in other tri- 
atomines as well as in many other blood-feeding arthro- 
pods including flies, mosquitoes and ticks (Mende et al. 
1999, Campos et al. 2002, Takac et al. 2006, Zhou et al. 
2006, Araujo et al. 2007, Mulenga et al. 2007, Ribeiro et 
al. 2007, Meiser et al. 2010). Kazal-type inhibitors are 
known to inhibit a range of serine proteinases. Native 
Kazals from blood-feeding arthropods inhibit thrombin, 
trypsin, factor Xlla, subtilisin A, elastase, chymotrypsin 
and plasmin (Friedrich et al. 1993, Campos et al. 2002, 
2004, Lovato et al. 2006, Meiser et al. 2010). 
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Kazal-type domains are characteristically 40-60 
amino acids long and inhibitors may contain single or 
multiple active domains. Six cysteine residues forming 
three disulfide bridges, C,:C 5 , C,:C 4 , C 3 :C 6 distinguish 
the conserved structure within classical and non-clas- 
sical Kazal-type domains. The predicted reactive site, 
PI amino acid residue, is located at position C,-X-P1 
and determines specificity within Kazal-type inhibitors 
(Kanost 1999). Within the domain, outside of the con- 
served cysteine residues, there are high amounts of vari- 
ability in other amino acid residues (Rimphanitchayakit 
& Tassanakajon 2010). 

Phlebotomine sandflies (Diptera: Psychodidae) are 
vectors of viruses, bacteria and parasites of the genus 
Leishmania. Transmission of Leishmania to suitable 
vertebrate hosts generally occurs during blood-feeding 
through the bite site of an infected sandfly vector [re- 
viewed by Ramalho-Ortigao et al. (2010)]. 

Midgut transcriptome analyses of Phlebotomus pa- 
patasi, the principal vector of Leishmania major, revealed 
two Kazal-type serine proteinase inhibitors, PpKzll and 
PpKzl2 (Ramalho-Ortigao et al. 2007). These were the 
first Kazal-type serine proteinase inhibitors identified 
from sandflies. The mature PpKzll cDNA is 231 base 
pairs (bp) encoding a 77 amino acid protein containing 
a single Kazal-type domain (GenBank ID: EU045342). 
The mature PpKzl2 cDNA is 267 bp encoding an 89 
amino acid protein (GenBank ID: JX171681). PpKzll 
and PpKzl2 have only 28% identity and 42% similarity 
in amino acid sequences (Ramalho-Ortigao et al. 2007). 
Both PpKzll and PpKzl2 have predicted signal peptides, 
suggesting that they are secreted in the midgut. 

We are interested in the role of these proteins in P. 
papatasi as inhibitors of serine proteinases and their 
potential effects on blood digestion. We have analysed 
deduced sequences of the PpKzll and PpKzl2 for pre- 
dicted activity and similarity, evaluated the expression 
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of PpKzll and PpKzl2 in developmental stages, adult 
female midguts and whole adult males and conducted 
in vitro analysis of inhibition activity of a recombinant 
PpKzl2 protein. 

MATERIALS AND METHODS 

Sandflies - P. papatasi Israel strain was reared in the 
Biology of Disease Vectors laboratory at the Department 
of Entomology, Kansas State University. Flies were main- 
tained on 30% sucrose solution at 27°C and 70% humid- 
ity with 12 h light and dark cycles. For blood feeding, 
sandflies were allowed to feed approximately 30 min on 
a BALB/c mouse anesthetised with 3 mg ketamine (Ket- 
aset, Fort Dodge Animal Health, Fort Dodge, IA, USA) 
and 0.12 mg xylazine (AnaSed, Acorn Inc, Decatur, IL, 
USA) per mouse (100 mg/kg of ketamine and 4 mg/kg of 
xylazine). Use of animals was preapproved by the Kan- 
sas State University Institutional Animal Care and Use 
Committee under protocols 2747, 2748 and 2749. Infec- 
tious blood meals contained L. major amastigotes and 
were offered artificially, while simultaneously a control 
set of sandflies were fed on uninfected blood as previ- 
ously described (Coutinho-Abreu et al. 2010a). 

At 20 h post-blood meal (PBM) all blood-fed flies 
were briefly anesthetised with C0 2 and examined under 
a dissecting microscope. Fully fed flies (i.e., abdomen 
fully distended) of similar size were selected for dis- 
section. Midguts were dissected in 30 |xL IX phosphate 
buffered saline RNAse free with ELIMINase (Fisher, 
Scientific, Pittsburgh, PA, USA) treated tools and equip- 
ment. Dissected midguts were then transferred to 50 
|xL of RNA later (Qiagen, Valencia, CA, USA), homog- 
enised with a hand-held homogeniser for approximately 
20 s and placed at -80°C. 

Sequence analysis - PpKzll and PpKzl2 were previ- 
ously identified from P. papatasi cDNA midgut libraries 
(Ramalho-Ortigao et al. 2007). Molecular weights and 
isoelectric points (pi) were predicted using the Swiss In- 
stitute of Bioinformatics ExPASy tools (Gasteiger et al. 
2003). Sequences similar to PpKzll and PpKzl2 were 
identified in National Center for Biotechnology Informa- 
tion using BLASTP for the non-redundant protein data- 
base (Altschul et al. 1997). The conserved six cysteine 



domain in PpKzll and PpKzl2 was used for multiple se- 
quence alignments (MSA) with selected sequences from 
blast results. Protein sequence alignments were performed 
using CLUSTALW2 (Larkin et al. 2007) and manual ed- 
its were performed in Jalview version 2 (Waterhouse et 
al. 2009). A Lutzomyia longipalpis Kazal2 contig (69116) 
was identified using BLAST searching for homologs of 
PpKzl2 in the L. longipalpis Lion 0.1 preliminary Ge- 
nome Assembly on the Baylor College of Medicine Hu- 
man Genome Sequencing Center website (hgsc.bcm.tmc. 
edu/project-species-i-Lutzomyia_longipalpis.hgsc). The 
sequence was translated with Swiss Institute of Bioinfor- 
matics ExPASy (Gasteiger et al. 2003). 

RNA extraction and cDNA synthesis - Total RNA 
was extracted from whole sample pools or individual 
dissected midguts using the RNeasy Mini Kit (Qiagen) 
and eluted in 40 |xL of RNase-free water. Three RNAs 
were obtained for each developmental stage from pools 
of 20 eggs, 10 Lj larvae and five each for stages L 2 , L y 
L 4 and pupae. Extracted RNA was treated with TURBO 
DNase (Ambion, Austin, TX, USA) to eliminate any re- 
sidual genomic DNA. Up to 100 ng of each RNA was 
used for first strand cDNA synthesis and was added to 
3.3 |xM oligo-dT 20 primer, 0.67 mM deoxynucleotide tri- 
phosphates and RNase-free water to total volume of 15 
uL. Samples were incubated at 65°C for 5 min and then 
placed on ice for 1 min. Addition of 4 |iL of 5X Super- 
Script III Reverse Transcriptase First-Strand Buffer, 5 
mM DTT, 0.5 |xL RNaseOUT (40 units/uL) and 1 |xL 
of Superscript III Reverse Transcriptase (200 units/|xL) 
(Invitrogen, Carlsbad, CA, USA) was followed with 1 h 
incubation at 50°C. All cDNA was stored at -20°C. 

Real-time polymerase chain reaction (RT-PCR) - 
PpKzll and PpKzl2 relative expression was analysed 
in non-blood-fed and blood-fed adult female sandflies. 
Individual midguts were dissected from non-blood-fed 
flies (0 h) and blood-fed flies at 24 h, 48 h and 72 h PBM. 
Total RNA was extracted from individual midguts and 
used for first-strand cDNA synthesis. RT-PCR was car- 
ried out on an Eppendorf Mastercycler ep Realplex 4 in 
8 uL reactions. Forward and reverse 0.3 |xM primers 
(Table) were mixed with 4 |xL iQ SYBR green Super- 



TABLE 
Complete list of primers 

Annealing 



Primer 


Primer sequence 5 '-3' forward 


Primer sequence 5 '-3' reverse 


(°C) 


PCR 


PpKzl859 


GCACCAGCCCAAAAGACC 


TCACTGCAATCTGATGGCGC 


56.5 


PCR 


VR1020 


ACAGGAGTCCAGGGCTGGAGAGAA 


AGTGGCACCTTCCAGGGTCAAGGA 


49 


PCR 


PpKzl2-R-His 


GCACCAGCCCAAAAGACC 


His tag'-CTGCAATCTGATGGCGC 


60 


PCR" 


PpKzll_137 


AGAGCGTTACCTGTCCTTG 


CCAGCGAATACTGAGGTTC 


58 


RT-PCR 


PpKzl2_152 


AATGAATGTCTGAAGGCCTG 


CCTTGGGATTTCACCTCCC 


58 


RT-PCR 


Pp40S_S3_136 


GGACAGAAATCATCATCATG 


CCTTTTCAGCGTACAGCTC 


58 


RT-PCR 



a: His tag-TCAGTGGTGATGGTGATGATG; b: touchdown polymerase chain reaction (PCR); RT: real-time. 
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mix (BioRad, Hercules, CA, USA) and added to 0.2 uL 
cDNA and 3.32 |iL molecular grade water (Invitrogen). 
All cDNA samples were run in duplicate for PpKzll 
and PpKzl2 and in parallel for 40S ribosomal protein S3 
(GenBank accession FG1 13203). Reactions were carried 
out 40 cycles of 95°C/30 s, 58°C/1 min and 72°C/30 s, 
followed by 95°C/15 s, 60°C/15 s and a melt curve up to 
95°C/20 min. C T values from the Realplex Software were 
used for expression analysis. 

Expression levels of mRNA were calculated with the 
comparative C T method as previously described (Coutin- 
ho-Abreu et al. 2010b). Briefly, C T values were normalised 
to the expression of a non-regulated internal control gene, 
40S ribosomal protein S3 and then normalised to a cali- 
brator. Calibrators for analysis of temporal, developmental 
and infected expression were mean averages of expression 
in 0 h, eggs and non-infected blood-fed samples respec- 
tively. Comparative C T method: AAC T = [AC T Variant X 
Sample] - [average (AC T Calibrator Samples)], where vari- 
ant X equals time points or tissue type. Fold change was 
calculated by 2 aAC c (Livak & Schmittgen 2001). Mean fold 
change of at least five individual samples or three pools 
were graphed for each time point or tissue. Distribution 
of the data was tested with the Kolmogorov-Smirnov test 
for normality and Levene's test for equality of variance. 
Nonparametric data was logarithmically transformed 
for statistical analysis. Data was evaluated with one-way 
analysis of variance and a parametric t test with the Bon- 
ferroni correction for multiple comparisons. For temporal 
expression profiles of L. major infected sandflies, statisti- 
cal analysis used two-tailed unpaired t tests for paramet- 
ric analysis and the two-tailed Mann-Whitney U test for 
nonparametric statistical comparisons. Prism 5 Software 
(GraphPad, La Jolla, CA, USA) was used for all graphing 
and statistical analysis. 

Recombinant protein expression and purification - 
The mature (minus signal peptide) PpKzl2 cDNA was 
amplified using the forward primer PpKzl859 and the 
reverse primer PpKzl2-R-His containing a 6X-His tag 
on its 3' end (Table), touchdown reverse transcriptase 
PCR was performed as follows, 95°C/3 min, three cy- 
cles of 94°C/1 min, 72°C/1 min, three cycles of 94°C/1 
min, 68°C/1 min, 72°C/1 min, five cycles of 94°C/1 min, 
62°C/1 min, 72°C/1 min, 25 cycles of 94°C/1 min, 60°C 1 
min, 72°C 1 min, finished with 72°C 5 min. Two microli- 
tres of the PCR product was separated on an agarose gel 
for analysis and to assess concentration. The mature Pp- 
Kzll was cloned into VR1020-TOPO vector as described 
previously (Ramalho-Ortigao et al. 2005, Oliveira et al. 
2006). Insert-containing clones were screened by PCR 
(Table) and orientation was confirmed by sequencing. 
Plasmid purification was as described by Oliveira et al. 
(2006). Final concentration was 2.5 mg/mL and plasmid 
sequence was confirmed by sequencing. 

The recombinant rPpKzl2 was expressed in CHO-S 
free style cells, following transfection using 37.5 ng of 
purified plasmid following the manufacturer's protocol 
(Invitrogen). Transfected CHO supernatant was collected 
after 72 h of culture, concentrated using a 3 kDa cut-off 
Centricon filter (Milipore, Billerica, MA, USA) and pu- 
rified by nickel-nitrilotriacetic acid chromatography with 



a gravity flow column. The column was washed with 15 
mL of 20 mM sodium phosphate buffer-300 mM sodium 
chloride-20 mM imidazole, eluted with 5 mL 20 mM so- 
dium phosphate buffer-300 mM sodium chloride-300 mM 
imidazole and the eluted rPpKzl2 was concentrated to 1.5 
(.ig/(.iL. Two hundred and fifty nanograms of protein were 
analysed by sodium dodecyl sulfate polyacrylamide gel 
electrophoresis using 4-12% reducing Bis-Tris NuPAGE 
pre-cast gel purchased from Invitrogen. The protein was 
transferred to nitrocellulose and incubated with anti-His 
antibody (Santa Cruz, Santa Cruz, CA, USA) overnight 
at 4°C and followed by three washes of 10 min each in tris 
buffered saline buffer with 0.1% Tween-20 (TBS-T). The 
blot was incubated with anti-mouse antibody conjugated 
to alkaline phosphatase (Promega, Madison, WI, USA) 
diluted 1:10,000 in TBS-T for 1 h at room temperature 
and washed in TBS-T. The protein bands were visualised 
using the Western Blue substrate (Promega). 

Inhibition assays - The inhibition activity of rPpKzl2 
was tested against human a-thrombin and trypsin and bo- 
vine a-chymotrypsin. Increasing concentrations of rPpK- 
zl2 were pre-incubated with 0.05 \iM human a-thrombin 
(Calbiochem, EMD Chemicals Inc, Gibbstown, NJ, 
USA), 2 uM trypsin (Sigma, St. Louis, MO, USA) or 0.25 
|j.M a-chymotrypsin (Calbiochem, EMD Chemicals Inc) 
in 50 mM Hepes-0.5% BSA, pH 7.3 for thrombin and in 
50 mM Tris-HCl, pH 8.0 for trypsin and a-chymotrypsin. 
Each enzyme and rPpKzl2 combination was incubated 
for 15 min at 37°C in a 96-well non-binding microtitre 
plate. Chromogenic peptide substrate H-D-Phenylalanyl- 
L-pipecolyl-Larginine-p-nitroaniline dihydrochloride 
(S-2238) (Chromogenix, diaPharma, West Chester Town- 
ship, OH, USA), Na-Benzoyl-D,L-arginine 4-nitroanilide 
hydrochloride (BAPNA) (Sigma) or N-Succinyl-L-ala- 
nyl-L-alanyl-L-prolyl-L-phenylalanine 4-nitroanilide 
(Suc-AAPF-pNA) (Sigma) was added at increasing con- 
centrations for a-thrombin, trypsin or a-chymotrypsin re- 
spectively for a total reaction volume of 100 \lL. Inhibiton 
of trypsin activity was measured for 3 nM, 30 nM and 
300 nM rPpKzl2 at increasing concentrations of BAPNA 
(25 nM, 125 uM, 250 yiM, 500 uM and 1000 jiM). Inhibi- 
tion of a-chymotrypsin activity was measured for 0.0005 
nM, 0.005 nM and 0.05 nM rPpKzl2 and inhibition of 
a-thrombin was measured at 0.5 nM, 3 nM and 300 nM 
rPpKzl2 at increasing concentrations 250 uM, 500 (0.M 
and 1000 fxM of Suc-AAPF-pNA or S-2238, respectively. 
The rate of proteinase hydrolysis of the chromogenic sub- 
strate was measured at 405 nm every 35 s during the reac- 
tion with a Biotek Synergy HT microplate reader (Biotek, 
Winooski, VT, USA). Each reaction was run in triplicate 
and each assay was repeated at least twice. 

Graphs of initial velocity (V) vs. substrate concentra- 
tion [S] were fit with the Michaelis-Menten equation to 
obtain the kinetic constant (K ) and maximum velocity 
(V ), v = Vmax[S] (Copeland 2000). Residual activity 
Km + [S] 

in the presence of different concentrations of rPpKzl2 
was calculated with apparent V values, residual activ- 

1 1 max 7 

ity = Vmax x 100 (Copeland 2005). 
V max, 0 
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RESULTS 

Sequence analysis - Both PpKzll and PpKzl2 code 
for six cysteine residues in a conserved arrangement 
characterised as a non-classical Kazal-type domain. 
Predicted molecular weights and isoelectric points for 
PpKzlland PpKzl2 are estimated to be 6.4 kDa and 5.22 
pi and 7.6 kDa and 6.10 pi respectively. In PpKzll an 
arginine residue is in the deduced PI site, the predicted 
active site for Kazal-type inhibitors (Fig. 1A) and PpK- 
zl2 contains a tyrosine in the PI site (Fig. IB). Arginine 
has been shown to confer thrombin and trypsin inhibi- 
tory activities and tyrosine in the PI commonly shows 
chymotrypsin inhibitory activity (Kanost 1999). 

PpKzll has 81% similarity and 73% identity to a puta- 
tive protein identified in the New World sandfly L. lon- 
gipalpis, vector of Leishmania infantum chagasi (Jochim 
et al. 2008, Pitaluga et al. 2009) (Fig. 1A). PpKzll has 
conserved sequence features previously described in non- 
classical Kazal-type domains in blood-feeding and non- 
blooding insects such as P-X-C 3 -G-X 4 -T-Y-X-N-X-C 4 and 
G-X-C 6 , with (X) representing various residues (Augustin 
et al. 2009). A MSA with the top blast results for PpKzll 
was assembled as described in Materials and Methods sec- 
tion and displayed high conservation of arginine in the PI 
site for this group of Kazals (Supplementary data). PpKzl2 
is also similar to another predicted protein in L. longipal- 
pis (Ramalho-Ortigao et al. 2007), but to a lesser degree 
with only 44% identity and 53% similarity (Fig. IB). While 
the amino acids in the predicted PI site in the P. papatasi 
and L. longipalpis proteins differ, tyrosine and phenylala- 
nine do share similar structural and chemical properties 
and have both been shown to inhibit chymotrypsin. The 
Kazal-type domains in the PpKzl2 MSA displayed large 
diversity in PI residues (Supplementary data). 

Non-classical Kazal-type domain patterns are par- 
tially conserved in PpKzl2 including regions P-X-C 3 and 
G-X-C 6 , (Fig. IB, Supplementary data). PpKzl2 also has 
more residues between C 3 and C 4 shifting the location of 
the fifth cysteine closer to the C-terminus, which has been 



seen in other non-classical Kazal-type domains (Hemmi 
et al. 2005, Rimphanitchayakit & Tassanakajon 2010). 
Conserved residues specific to the PpKzl2 MSA include 
N-C -E/Q and a phenylalanine located four residues up- 
stream of the fourth cysteine (Supplementary data). 

Expression profiles - Expression of PpKzll and Pp- 
Kzl2 in the female midgut increased after blood feeding. 
Temporal expression was analysed 0 h, 24 h, 48 h and 72 
h PBM. PpKzll transcript expression was up-regulated 
at 24 h and 48 h PBM (p < 0.05, p < 0.001) (Fig. 2A). 
After a significant increase in expression at 48 h PBM, 
PpKzll expression decreased to pre -blood feeding levels 
(0 h) between 48-72 h PBM (p < 0.01). Expression of 
PpKzl2 was up-regulated 24 h, 48 h and 72 h PBM (p 
< 0.01, p < 0.001, p < 0.05) (Fig. 2B). Transcript levels 
were up-regulated at 24 h and continued to increase sig- 
nificantly at 48 h PBM (p < 0.01). PpKzl2 expression was 
then down-regulated by 72 h (p < 0.001) with expression 
at 72 h decreasing to levels similar to 24 h expression. 

Following results indicating that PpKzll and PpK- 
zl2 expression is regulated following a blood meal, we 
then investigated if these transcripts are expressed in 
developmental (non-blood feeding) stages. Expression 
profiles of developmental stages for PpKzll and PpKzl2 
show both transcripts expressed during early develop- 
ment (Fig. 3). Both PpKzll and PpKzl2 are expressed 
in larval stages h v L 2 , L 3 , L 4 and pupa at constant levels 
showing no significant differential regulation in expres- 
sion between developmental stages. Both PpKzll and 
PpKzl2 are expressed in whole male tissues, but expres- 
sion was not detected in eggs for either Kazal transcript 
(Supplementary data). 

PpKzll and PpKzl2 expression was further analysed 
at 24 h, 48 h and 72 h following an infective blood feed- 
ing with 5 x 10 6 L. major amastigotes per mL of blood. 
No significant difference in the mRNA expression lev- 
els of PpKzll and PpKzl2 between non-infected vs. L. 
major infected flies were detected in these three time 
points (Fig. 4). 
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Fig. 1: PpKzll and PpKzl2 alignments with Lutzomyia longipalpis Kazal domains. PpKzll (GenBank ID: EU045342) (A) and PpKzl2 (Gen- 
Bank ID: JX171681) (B) are both similar to putative proteins in L. longipalpis with Kazal-type domains: LIKzll (GenBank ID: ABV60319) and 
LlKzl2 (contig 69116). Conserved residues are in black and similar residues are in grey. Predicted signal peptides are underlined, asterisks mark 
predicted PI residues, conserved cysteines are marked (C) and gaps are indicated by dashes. 
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Inhibition assays - Inhibition activity of rPpKzl2 was 
tested for a-thrombin, trypsin and a-chymotrypsin en- 
zymes. Residual activity of enzymes in the presence of 
rPpKzl2 was reduced to 9.4% for a-chymotrypsin, 33.5% 
for a-thrombin and 63.9% for trypsin (Fig. 5). Both V 
and K decreased in all inhibition assays with increas- 
ing concentrations of rPpKzl2 (Supplementary data). 
Recombinant PpKzl2 inhibited a-chymotrypsin at the 
nanomolar level and inhibited a-thrombin and trypsin at 
micromolar levels. 
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Fig. 2: PpKzll and PpKzll expression in adult females post-blood 
meal (PBM). PpKzll and PpKzl2 mRNA expression levels are regu- 
lated after a blood meal. A: PpKzll is up-regulated 24 h and 48 h 
PBM with highest expression at 48 h PBM. By 72 h expression is 
down-regulated to levels similar to 0 h; B: PpKzll is up-regulated 24 
h, 48 h and 72 h PBM. Expression is highest at 48 h and decreases be- 
tween 48-72 h PBM. Values are the mean fold change of five or more 
individual midguts with standard error of the mean. Expression was 
calibrated to 0 h expression levels. Analysis used ANOVA t test with 
the Bonferroni correction for multi-comparisons. *: p < 0.05; **: p < 
0.01; ***: p < 0.001. 



DISCUSSION 

Kazal-type inhibitors are a diverse group of serine 
proteinase inhibitors with a wide range of roles in in- 
vertebrates. In blood-feeding triatomines, Kazal-type 
inhibitors in the midgut prevent coagulation of the blood 
meal (Friedrich et al. 1993, Mende et al. 1999, Campos et 
al. 2002, 2004, Araujo et al. 2007, Meiser et al. 2010). 

Here, we characterised two single domain non-clas- 
sical Kazal-type inhibitors from the sandfly P. papatasi. 
PpKzll and PpKzl2 mRNA transcripts are expressed in 
non-blood-fed and blood-fed female midguts and expres- 
sion is regulated by the blood meal with up-regulation at 
24 h and 48 h PBM. The decrease in PpKzll and PpKzU 
expression detected around 72 h PBM correlates with 
the completion of blood meal digestion, which culmi- 
nates with the midgut emptying between 72-144 h PBM. 
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Fig. 3: PpKzll and PpKzll expression in larval stages and pupa. A: 
PpKzll was expressed in all larval stages and pupae. PpKzll expres- 
sion was not significantly different when compared between larval 
stages; B: PpKzll was also expressed in all larval stages and pupae 
at similar expression levels. Five or more individuals were pooled for 
each developmental stage and this was repeated for a total of three 
replicates. Values are the mean fold change with standard error of the 
mean. Expression was calibrated to expression levels in eggs. ANO- 
VA t test with the Bonferroni correction for multi-comparisons was 
used for statistical analysis. L: larval stage; P: pupa. 
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Furthermore, the expression levels of both PpKzIl and 
PpKzl2 remain constant between 72-144 h PBM (Supple- 
mentary data). Such expression profiles of PpKzll and 
PpKzl2 are suggestive of a role in digestion for their re- 
spective proteins. In addition, as PpKzll and PpKzl2 also 
are expressed in all larval stages, pupae and males, inhi- 
bition during digestion is likely not specific to serine pro- 
teinases involved in the coagulation cascade, but rather 
serine proteinases engaged across life stages and sexes. 

The predicted PpKzll is similar to a single domain 
non-classical Kazal-type inhibitor from Aedes aegypti, 
AaTI (Ribeiro et al. 2007). Interestingly, a recombinant 
AaTI was shown to inhibit trypsin and plasmin, with 
weak inhibition of thrombin activity; the AaTI transcript 
also was shown to be expressed in larva, pupa, male and 
female tissues (Watanabe et al. 2010, 2011). PpKzll is 
also similar to the multi-domain Kazal-type inhibitors 
infestin and dipetalogastin, identified in T. infestans 
(Campos et al. 2002), and Dipetalogaster maximus 
(Mende et al. 2004), respectively, but with the highest 



identity to infestin's domain-4. This domain was found 
to strongly inhibit factor Xlla, plasmin and trypsin, with 
no activity for thrombin (Campos et al. 2002, 2004). 
Consistent with previous findings, PpKzll as a non-clas- 
sical Kazal-type domain displays a predicted active site 
residue that suggests it likely possess inhibitory activity 
for trypsin-like serine proteinases. 



Chymotrypsin 




O 



PpKzll 



11 <W> 

7"T OOO 



24 h 



I Nl 

48 h 



I Nl 

72 h 




rPpKztt nM 



Trypsin 



rPpKzl2 nM 



PpKzl2 





AA* «S» 




A**A O 

A O 


Tp- <>oo 

A 


A 


o 

A 






1 1 
1 Nl 


1 1 
1 Nl 


1 1 
1 Nl 


24 h 


48 h 


72 h 



Fig. 4: PpKzll and PpKzll expression in adult females infected with 
Leishmania major. Temporal expression profiles 24 h, 48 h and 72 h 
post-infective blood meal (I) (A) and post-non-infected blood meal 
(Nl) (O). Eight individual midguts were assayed for each infected 
and non-infected time point. PpKzll and PpKzU expression was not 
significantly different 24 h, 48 h and 72 h I when compared to Nl 
control groups. Bars are the mean fold change of eight individual 
midguts. Expression was calibrated to expression in Nl controls. Sta- 
tistical analysis used two-tailed unpaired t tests and two-tailed Mann- 
Whitney U tests for parametric and nonparametric comparisons re- 
spectively (p < 0.05). 




Thrombin 



rPpKzl2 nM 

Fig. 5: rPpKzl2 enzyme inhibition activity. Activity was measured at 
increasing concentrations of both rPpKzl2 and substrate. Reactions 
were fit with Michaelis-Menten non-linear regression and apparent 
maximum velocity (V ) values were used to calculate residual activ- 
ity. Inhibition of a-chymotrypsin activity was observed with decreas- 
ing V . Activity of 0.25 uM a-chymotrypsin was reduced to 9.4%. 
Residual activity of 2 uM trypsin was reduced to 63.9%. Activity 
of 0.05 u.M a-thrombin in the presence of rPpKzl2 was reduced to 
33.5%. Reactions were run in triplicate and each graph represents one 
of two replicates of each experiment. 
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PpKzl2 on the other hand is similar to Kazal-type 
domains from dipteran, lepidopteran and hymenopteran 
species. Though no functional characterisation for these 
Kazal domains have been described, putative proteins 
were identified in expressed sequence tag and cDNA li- 
braries of immune-challenged insects (Bartholomay et 
al. 2004, Gandhe et al. 2006). 

A recombinant PpKzl2 was obtainedand tested against 
various substrates. Inhibition activity of rPpKzl2 was ob- 
served for a-chymotrypsin, a-thrombin and trypsin, in 
agreement with previous reports on single-domain Ka- 
zal-type inhibitors having activity against multiple serine 
proteinases (Nirmala et al. 2001, Watanabe et al. 2010). 
The ability of PpKzl2 to inhibit serine proteinases in P. 
papatasi midgut is dependent upon the rate of inhibition 
and concentrations present in the midgut (Kanost & Jiang 
1996) and therefore in vivo activity may be enzyme spe- 
cific. Whereas rPpKzl2 inhibited a-thrombin, the inhibi- 
tion activity for a-chymotrypsin was the strongest. We 
previously characterised two chymotrypsin-like and four 
trypsin-like proteases from P. papatasi and demonstrat- 
ed that chymotrypsin and trypsin activities in the midgut 
of this sandfly peak between 27-48 h PBM and by 72 h 
PBM no such activities were detected (Ramalho-Ortigao 
et al. 2003). Also, as our results indicate, the peak in RNA 
abundance for Kazals in P. papatasi is 48 h PBM. These 
data, together with the observations that rPpKzl2 inhib- 
ited both chymotrypsin and trypsin and expression of the 
mRNA was also observed in non-blood-feeding life stag- 
es, suggest to us that PpKzl2 is more likely involved in 
regulating digestive proteases than blood fluidity within 
the midgut. Knock down by injection of 127 ng of double 
stranded RNA produced against each target did not affect 
mRNA expression levels of PpKzU and PpKzll in the 
midgut of P. papatasi and therefore analysis of effects 
on blood meal digestion rate via haemoglobin levels in 
female midguts were not informative (data unpublished). 

Some Kazals have been shown to have immune-like 
activity; however there was no response in transcript ex- 
pression of PpKzll and PpKzU during L. major infection. 
No effects were observed on PpKzll and PpKzl2 expres- 
sion during L. major infection in the midgut at 24 h, 48 h 
or 72 h post-infective-blood meal. It has been described 
in sandflies that infection with Leishmania leads to mod- 
ulation of trypsin-like activity in the midgut during di- 
gestion, suggesting that modulation of trypsin activity al- 
lows the parasites to survive (Borovsky & Schlein 1987, 
Sant Anna et al. 2009, Telleria et al. 2010). This has been 
supported with data showing that RNAi of a trypsin gene 
increased parasite numbers during infection (SantAnna 
et al. 2009). The dynamics of serine proteinases and ser- 
ine proteinase inhibitors in the midgut are not only cru- 
cial to sandfly metabolism and digestion, but may also 
affect Leishmania development. Further characterisation 
of the serine proteinase cascades and their inhibitors in 
P. papatasi may provide insight into the complex interac- 
tions that constitute vector competence. 
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PpKzll multiple sequence alignment. Phlebotomus papatasi PpKzll (GenBank ID: EU045342), Culicoides sonorensis (GenBank ID: 
AAV84258), Drosophila yakubaa (GenBank ID: XP_002088400), Anopheles darlingia (GenBank ID: ACI30143), Lutzomyia longipalpis (Gen- 
Bank ID: ABV60319), Aedes aegypti AaTi (GenBank ID: ABF18209), Culex quinquefasciatus (GenBank ID: XP 001868221), Ochlerotatus 
triseriatus salivary (GenBank ID: ACU30983), Aedes albopictus (GenBank ID: AAV90671), Triatoma infestans infestin domain 4 (GenBank 
ID: AAK57342), Anopheles gambiaea (GenBank ID: XP 001230687), An. gambiae b (GenBank ID: XP 317819), An. gambiaea (GenBank 
ID: EAA12788), Drosophila yakubab (GenBank ID: XP 002088399), D. yakubac (GenBank ID: XP 002088401), D. yakuba i (GenBank ID: 
XP 002088397), Drosophila pseudoobscura pseudoobscuraa (GenBank ID: XP 001356962), D. pseudoobscura pseudoobscurajo (GenBank 
ID: XP 001356963) and An. darlingijo (GenBank ID: ACI30165). Asterisk means the predicted PI residue, gaps are indicated by dashes, con- 
served cysteines are in black and residues with more than 50% conserved identity are in shades of grey. 
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PpKzl2 multiple sequence alignment. Phlebotomus papatasi PpKzl2 (GenBank ID: JX171681), Anopheles darlingi (GenBank ID: ACI30205), 
Drosophila mojavensis (GenBank ID: XP 002000106), Culex quinquefasciatus (GenBank ID: XP 001842298), Aedes aegypti (GenBank ID: 
XP 001658905), Manduca sexta (GenBank ID: AAF16698), Nasonia vitripennis (GenBank ID: XP 001600330), Bombyx mori (GenBank ID: 
NP 001040250), Bombus terrestris (GenBank ID: XP 003401213), Drosophila pseudoobscura pseudoobscura (GenBank ID: XP 001359513), 
Drosophila persimilis (GenBank ID: XP 002017393), Tribolium castaneum (GenBank ID: XP 974370), Drosophila willistoni (GenBank ID: 
XP 002070657), Drosophila grimshawi (GenBank ID: XP 001994337), Solenopsis invicta (GenBank ID: ADC34234), Drosophila virilis (Gen- 
Bank ID: XP 002053264), Acromyrmex echinatior (GenBank ID: EGI69242), Harpegnathos saltator (GenBank ID: EFN89909), H. saltator 2 
(GenBank ID: EFN81812), Bombus impatiens (GenBank ID: XP 003486913), Apis florae (GenBank ID: XP 003692060), Camponotus f/ori- 
danus (GenBank ID: EFN62548) and Lutzomyia longipalpis (contig 69116). Asterisk means the predicted PI residue, gaps are indicated by 
dashes, conserved cysteines are in black and residues with more than 50% conserved identity are in shades of grey. 
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PpKzll and PpKzl2 expression in Phlebotomus papatasi males and not in eggs. Reverse-transcriptase polymerase chain reaction (PCR) was 
carried out in 20 pL reactions with cDNA from one whole male and a pool of 10 eggs on an Eppendorf Mastercycler gradient. Reactions were 
prepared with 10 pL of 2X GoTaq PCR master mix (Promega, Madison, WI, USA), 0.2 pM forward and reverse primers, 1 pL cDNA and mo- 
lecular grade water (Invitrogen, Carlsbad, CA, USA) to a total volume of 20 pL. Reactions were carried out as follows: 95°C/1 min followed by 
26 cycles of 94°C/30 s, 56.5°C/30 s and 72°C/1 min and a final step at 72°C/5 min. The primers PpKzll 1 1 , PpKzl859 and PpTubl48 specific for P. 
papatasi VpKzIl, PpKzI2 and (3-tubulin were used for PCR. PCR fragments (10 pL) were separated on a 1.5% agarose gel stained with ethidium 
bromide. Primers for PpTubl48 forward: GCGATGACTCCTTCAACAC and reverse: GTGATCAATTGTTCGGGATG. 
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Michaelis-Menten non-linear regression of rPpKzl2 inhibition. Initial velocity (V) over substrate concentration (S) was fit with Michaelis- 
Menten non-linear regression for each concentration of rPpKzl2. A reduction in maximum velocity and kinetic constant values was observed 
with increasing rPpKzl2 when compared to the fit of 0 nM rPpKzl2. BAPNA: Na-Benzoyl-D,L-arginine 4-nitroanilide hydrochloride; S-2238: 
H-D-Phenylalanyl-L-pipecolyl-Larginine-p-nitroaniline dihydrochloride; Suc-AAPF-pNA: N-Succinyl-L-alanyl-L-alanyl-L-prolyl-L-pheny- 
lalanine 4-nitroanilide. 
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PpKzll and PpKzl2 expression 72-144 h post-blood meal (PBM). Temporal analysis with semi-quantitative reverse transcriptase-polymerase 
chain reaction (PCR) indicated that PpKzll and PpKzll transcript expression remains constant 72-144 h PBM. Time points were pools of five 
midguts from female sandflies. PCR was carried out in 20 uL reactions on an Eppendorf Mastercycler gradient. Reactions were prepared with 10 
uL of 2X GoTaq PCR master mix (Promega, Madison, WI, USA), 0.2 pM forward and reverse primers, 1 pL cDNA and molecular grade water 
(Invitrogen, Carlsbad, CA, USA) to a total volume of 20 pL. Reactions were carried out as follows: 95°C/1 min followed by 26 cycles of 94°C/30 
s, 56.5°C/30 s and 72°C/1 min and a final step at 72°C/5 min. The primers PpKzll 11, PpKzl859 and PpTubl48 specific for Phlebotomus papatasi 
PpKzll, PpKzll and (3-tubulin were used for PCR. PCR fragments (10 pL) were separated on a 1.5% agarose gel stained with ethidium bromide 
alongside 0.5 pg of exACTGene cloning DNA Ladder (Fisher, Scientific, Pittsburgh, PA, USA). Intensities of PCR fragments were standardised 
to (3-tubulin and compared with known quantities of the reference ladder using Total Lab 100 software (BioSystematica, Sarnau, UK). 



